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INTRODUCTION 


The  late-time,  fully-developed  stage  of  sausage-like  "beading"  of 
z-pinch  discharges  can  give  rise  to  enhanced  electron  heating  (and  thus 
enhanced  radiative  losses) ,  caused  by  Ohmic  anomalous  heating  in  the 
constricted  regions.  In  this  report,  this  transfer  of  energy  from  conden¬ 
sations  of  magnetic  field  energy  to  radiation  is  examined  quantitatively, 
based  on  a  simplified  dynamical  model  for  the  nonlinear  instability 
driving  terms. 


'*>  Constricted  portions  of  the  discharge,  with  low  density  and  cross- 
section,  and  mostly  or  entirely  "anomalous”  current,  alternate  with  the 
higher  density  "^eads"^  which  carry  current  classically.  The  extreme 
limit  of  this  phenomenon  is  that  of  multiple  diodes  in  series,  with  the 
nearly-evacuated  low  density  regions  considered  as  bipolar-flow  diodes, 
with  pinched  electron  flow,  as— suggested  by  Goldstein  (1.9 78£>  In  all 
probability,  the  low-density  regions  cannot  evacuate  to  the  extent  required 
for  such  vacuum-diode  behavior. 


^  >  The  overall  resistive  heating  rate  is  of  course  VI,  with  the  total 
current  I  given  by  appropriate  circuit  equations,  but  the  local  heating 
rates  for  electrons  in  the  low  density  regions  are  balanced  by  increased 
radiative  loss  when  these  hotter  electrons  collide  with  the  denser  blobs 
of  plasma.  The  blobs  cannot  respond  hydrodynamically  to  the  increased 
heating  before  radiation  loses  the  deposited  energy. 


CURRENT 


We  first  examine  the  conditions  under  which  all  the  current  in  the 

low-density  regions  would  be  drift-speed  limited,  i.e.,  (approximately) 

when  the  classical  drift  speed  Vj  *  aE/nge  would  exceed  the  sound  speed 

c  .  From  a,  *  w2/4rrv  ,  and 
s  -*•  p  ei 

v^Uec-1 4 *)  -  0.9  X  10”  ^  (fg)1 *  (g^r)  T-^2  (1) 

we  get 

- l-8  * 10"  (is)  (r)2  (^7)  Tk2v2  w/o  (« 

where  z  is  the  degree  of  ionization  (ng  »  zn^) ,  Oj_  is  the  axial  conduc¬ 
tivity,  (across  magnetic  field  lines)  and  ZnA  is  the  plasma  parameter. 

The  classical  drift  speed  rises  rapidly  with  T  as  the  plasma  becomes  less 
collisional. 


This  is  to  be  compared  with 


(3) 


For  10-times  ionized  aluminum  (A  =  26,  z  ■  10),  one  sees  that  with 
electrons  at  keV  temperatures  the  classical  drift  speed  greatly  exceeds 
the  sound  speed  for  MV/cm  fields  unless  n^  is  of  order  1021  cm  or 
greater.  Elsewhere  we  argue  that  the  actual  drift  speed  cannot  much 
exceed  the  sound  speed,  so  that  the  current  is  drift-speed  limited  to  a 
value  approximately 


Ia(MA) 


(4) 


in  the  anomalous  ("a")  low-density  regions  of  cross-sectional  area  irr2. 

Since  the  current  density  is  divergenceless ,  this  must  be  the  current 

in  the  circuit. 


VOLTAGE  DROP  IN  LOW  DENSITY  REGIONS 


We  now  envision  a  number  N  of  such  millimeter  size  low-density 

"anomalous"  regions  of  length  l  alternating  with  denser  "classical" 

regions  ("beads")  of  length  &c.  Each  "classical"  region  has,  in  general, 

a  core  which  carries  no  current  because  it  has  not  been  penetrated  by  the 

magnetic  field  (and  hence  has  no  inductive  E) ,  a  classical  skin  current 

layer,  and  an  "anomalous"  (drift-speed  limited)  outer  corona.  The 

electrical  load  is  thus  envisioned  as  a  transmission  line  consisting  of 

N  elements  in  series,  with  each  element  as  in  figure  1.  Since  the  total 

current  is  divergenceless,  it  is  given  by  Eq.  (4),  and  depends  on  the 

applied  voltage  only  through  the  temperature.  If  the  corona  carries  only 

a  small  fraction  (1-0«1)  of  the  current,  the  DC  voltage  drop  along  each 

classical  bead  is  R  I  where  R  is  the  classical  resistance  of  the  bead, 

c  c 


R 


c 


i.e.  R(fl) 
c 


3.5  x  1<T4 


£.c(mm) 

w(rJ  -  r^Kmm*) 


(5) 


with  r£  the  bead  radius  and  ri(<  rfi)  the  current  penetration  radius. 

•  • 

One  can  see  that  when  inductive  (LI  and  LI)  effects  are  small  (l*e.» 
if  ia  0  when  I  peaks  or  dips) ,  the  classical  resistance  of  a  few  tens  of 
1  mm  beads  at  roughly  1  keV  is  quite  small  if  the  current  penetrates  the 
beads,  and  almost  all  the  voltage  drop  occurs  in  the  anomalous  regions. 

On  the  other  hand,  r, /r  tends  to  be  frozen  in  at  early  times  by  the 

c  -1 
formation  of  the  high  conductivity,  and  tends  to  have  values  of  order  10 

or  less,  so  the  resistance  is  larger  than  the  current-penetrated  value  by 

an  order  of  magnitude  or  more.  Still,  this  gives  only  a  small  resistance 

for  the  "classical"  portion  of  the  circuit  and  most  of  the  DC  voltage  drop 

would  occur  in  the  anomalous  regions. 
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JOULE  HEATING 

Setting  aside  for  now  the  question  of  inductive  corrections,  this 
allows  us  to  estimate  the  Joule  heating  rate  in  the  anomalous  regions  in 
the  DC  case: 

V  T  VI  Vc  n  t2 

E‘I  *  N  l  ~  N  l  c1  C6) 

a  a  a  a 

where  the  last  term  is  neglected  and  where  I  is  given  by  Eq •  (4) .  If  all 
this  energy  went  into  electron  temperature  and  none  escaped  to  ions  or 
photons,  it  would  correspond  to  a  rise 

•  2  VI  1  ,,v 

3  N  Jl  Trr2  zn. 
a  a  a  i 


T(kev/M)  »  15  VOW)  [snfc^lf  Tof  T; 


h 

keV  ’ 


a  number  of  order  15  keV/ns. 

The  transit  time  of  a  typical  electron  from  one  bead  to  the  next  at 
a  constant  sound  speed  c3  (T^y)  would  be 


T  ■  l  /c  *4.5x10  sec  x  i  (mm) 
as  a 


(fen 


or  a  few  ns,  so  that  an  electron  could  acquire  several  tens  of  keV  before 
plowing  into  the  adjacent  dense  bead  and  losing  its  energy  to  dense- 
plasma  thermal ’ energy  and  not-so-cold-target  bremsstrahlung.  In  fact,  a 
certain  amount  of  this  energy  is  lost  to  ions  and  radiation  during  transit 
across  the  anomalous  region,  and  we  now  give  estimates  showing  that  loss  to 
be  dominated  by  classical  heat  conduction. 
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COOLING  RATES  OFFSETTING  JOULE  HEATING 


The  classical  collisional  electron  cooling  time  to  ions  of  atomic 
mass  A  (taken  from  Book,  1980)  is 


T(keV/ns)  -  -  3.6  *  Hf2  (fj)’  (f )  (j^t)  t£t 

which  is  negligible  at  1  keV  and  above  (unless  ng  >  1021  cm--3) . 
radiative  cooling  gives  approximately 


T(keV/ns)  ~ 


(10) 


The 


(11) 


for  aluminum  with  Te  >_  1  keV  (Terry  and  Guillory,  1980b,  p.  5  and  plasma 
radiation  transport  computations  done  by  Duston  (1980)). 


The  classical  heat  flux,  which  is  dominated  by  the  Hall  flux  term 
B  x  V.  T  since  v  /u)  is  small  except  very  near  the  axis,  would  give  a 

^  6  0  C6 

cooling  rate 


T(keV/ns) 


where  r^  is  the  radial  temperature  scaleheight. 


The  modified  two-stream  instability  (Krall  and  Liewer,  1971;  Ott 
et  al,  1972)  grows  at  a  rate  about  half  the  lower-hybrid  frequency 


luLH 


/w  u  ~  6  x  1011  (~)  (—?)(.- -a  )  T. 

ce  ci  \A  /  \10/  \10l8/\l  nun  / 


H 

keV 


(13) 


(based  on  Eq.  (4)  for  I)  and  thus  can  easily  saturate  in  less  than  1  ns, 
leading  to  alteration  of  the  diffusive  heat  flux  as  well  as  being  the 
source  of  the  drift-speed  limiting  assumed  at  the  outset. 
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But  the  alteration  of  the  heat  flux,  according  to  the  form  of 
anomalous  heat  conductivity  given  by  Manheimer  (1979) , 

<  t+>  cT/16  eB 

is  small  and  gives  a  contribution 

TOceV/na)  -  l<f3  (f)  W  (J)  ] 


to  the  electron  cooling  in  transit  across  the  anomalous  low-density  regions. 
(The  inverse  dependence  on  n^  comes  only  from  the  linear  dependence  of  I , 
and  thus  B,  on  n^.) 

The  instability  tries  to  stabilize  by  bringing  a  small  number  of 
ions  up  to  the  sound  speed  (Vekshtein  et  al,  1971)  but  this  is  an  important 
energy  drain  for  the  electrons  only  where  Tg  »  Ti>  [For  the  time  develop¬ 
ment  of  T.  due  to  this  anomalous  heating,  one  may  consult  Lampe  et  al, 

(1971)3-  A  fraction  f.  ^  (m  /c.)  of  the  ions  is  accelerated  to  a  speed 

l  e  i  ^ 

comparable  with  cg,  i.e. ,  to  energies  Tg.  (The  directions  are  isotropized 
by  the  magnetic  field.)  There  is  thus  a  net  loss  of  ion  thermal  energy 
due  to  the  exit  of  these  particles  from  the  low  density  region,  which  is 
only  partly  compensated  by  the  thermal  influx  of  ions  at  temperature  T^. 

The  net  heat  loss  flux  per  cm2  of  area  is 


*l£i[(is)(lIe)-(i’'Thi)(!Ii)] 


(v^^  is  the  ion  thermal  speed);  and  when  the  instability  is  saturated 
and  quasi-steady  this  energy  loss  is  made  up  by  energy  transfer  from 
electrons  to  ions  via  unstable  waves: 


2  n.*e  (|  nl£i)[cs£e  -  vThl  Ti] 
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This  gives  rise  to  an  electron  cooling  rate 


ie0»v/»s)  <  .5  *  10-3  (ff  (ip)  x3«  .  (17) 

This  rate  is  small  because  of  the  small  number  of  ions  involved, 

>  relative  to  the  number  of  electrons.  In  fact,  for  nominal  parameters 

all  of  the  electron  cooling  rates  CIO)  -  (17)  are  small  compared  with 
the  Ohmic  heating  rate,  and  so  most  of  the  Ohmic  energy  gain  is  delivered 
to  the  blobs  after  transit  across  the  anomalous  region. 

I 


» 


» 


» 
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POWER  DELIVERY 


Assuming  the  dense  blobs  to  the  thicker  than  one  range  for  the  15-50 
keV  electrons  accelerated  into  them,  the  total  net  power  delivery  to 
them  is 


I 


» 


I 


» 


P«  4  (I  -  T  )  N7rr2  n  c 
2  e  eo  a  e  s 


(18) 


where  Tg  -  Teo  is  the  temperature  gain  by  Ohmic  heating,  less  losses, 
during  transit  from  blob  to  blob.  If  the  cooling  during  transit  is 
negligible  and  if  almost  all  the  voltage  drop  occurs  across  the  low-density 
regions,  then  the  power  delivered  is 


PCW)  -  0.36  (^)  (^*r)  (|5)3/2  (ff  (j£r) 


keV 


(19) 


with  n^,  T^eV  representing  either  their  average  values  of  their  initial 
(unheated)  values  in  the  low-density  region.  If  less  than  the  full  voltage 
drop  occurs  across  the  anomalous  regions,  the  value  of  V  above  is  adjusted 
downward  accordingly. 


I 


» 
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PEAK  ELECTRON  TEMPERATURE 


1 


I 


I 


» 


I 


I 


t 

s 

f 

r 
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In  the  low-density  regions,  the  electron  temperature  varies 
axially,  increasing  in  the  direction  of  electron  flow  because  of 
collisional  and  unstable-wave  isotropization  of  the  energy  gained 
from  the  axial  electric  field.  The  temperature  begins  at  approxi¬ 
mately  its  value  in  the  adjacent  classical  bead  (a  quasi-balance  of 
deposition  and  radiative  loss) ,  increases  as  the  electron  traverses 
the  low  density  region,  and  reaches  its  maximum  just  before  deposition 
i  the  next  "bead"  downstream. 

Since  the  heating  during  transit  has  the  approximate  form  T  « 
(from  Eq.  (8)),  one  has  for  the  case  of  negligible  heat  conduction  and 
radiation 


c8(t) »  ca(0)  +  cgt 

z(t)  -  c  (t) 
s 


(20) 


So  that  z^t)  =  c  (0)t  +  t  c  t2 
S  l  s 


and 


the  final  temperature  is  such  that 


c  (t.)  =  TP  +  2  c  l 

st  so  s  a 


(21) 


i.e.  Te(£a)  -  Teo  +  67°  keV  j^V(MV)  -  Llj  (22) 

from  Eqs.  (3)  and  (8),  as  an  overestimate  of  the  peak  temperature. 

(Na,  the  number  of  anomalous  regions,  is  observed  to  be  typically 
of  the  ordef  of  20  to  30,  while  the  inductively  corrected  voltage 
V  -  LI  is  one  or  two  MV.)  The  voltage  drop  along  the  classical  beads 
has  been  neglected  compared  with  that  over  the  anomalous  regions 
between  beads,  and  including  it  reduces  the  effective  voltage  appearing 
in  Eq.  (22)  to  that  occurring  only  over  the  N&  anomalous  regions.  In 
this  approximation  the  temperature  increases  linearly  with  distance 
across  each  acceleration  region. 
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The  indication  that  the  peak  electron  temperature  may  be  tens 
of  keV  implies  that  Bremsstrahlung  and  inner  shell  x-rays  can  be  produced, 
predominantly  on  the  cathode  ends  of  the  classical  beads  if  their  density 
is  high  enough,  or  throughout  the  beads  if  they  are  only  about  one  range 
thick  for  the  energetic  electrons. 
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INDUCTIVE  EFFECTS 


The  generalized  Ohm's  law  in  a  classical  medium  takes  the  form 

E»nJ  +  kj--xB  (23) 

e 

or  V  =*  Rl  +  LI  +  LI  .  (24) 


(k  is  a  constant  related  to  the  electron  inertia  and  the  geometry,  and 
L  is  the  self -inductance- )  Voltages  induced  by  changing  the  current- 
carrying  radius  r^  are  of  order 


E (MV/cm) 


21  (MA)  —■  (ns_1)  , 
a 


(25) 


which  becomes  of  .rder  1  MV  just  before  the  peak  compression.  The 
corresponding  impedance 


L  (Ohms)  2£(cm)  r  /r  (ns-1) 

a  a 


(26) 


has  maximum  value  on  the  order  of  1  Ohm,  in  series  with  1-4  Ohms  effec¬ 
tive  resistance,  and  so  is  not  a  negligible  correction  except  for  a  few- 

nanosecond  intervals  when  r  «0,  In  SPLAT  code  runs,  the  L  term  is 

a 

even  somewhat  larger,  at  its  maximum,  than  the  Ohmic  resistance  of  the 
classical  plasma  core.  But  where  the  current  becomes  limited  by  anomalous 
regions  in  series  with  classical  ones,  the  effective  resistance  R*  in 
the  equation 


V  *  R*I  +  .LI  +  LI 


(27) 


is  larger  (and  voltage-dependent,  as  V);  so  the  L  terra  no  longer 
dominates,  although  it  is  not  strictly  negligible  except  very  near  peak 
compression  and  for  the  "pause"  phase  of  "pause-and-collapse"  compressions. 
For  our  simple  description  we  neglect  the  L  term,  which  is  acceptable 
when  V/I  is  much  larger  than  L,  i.e., 
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(ns*1) 


«  1.4 


V(MV) 

Jl(cm) 


/A  \  ^  /10\3/2  ao‘*\  -1/2  {lma2\  .... 

(w  lr)  {—)\ev  <28> 


The  LI  term  may  be  governed  by  nonlinear  sausage  instability 

dynamics  not  yet  well  understood;  when  I  is  given  by  Eq.  (4),  I  is 
d  2 

proportional  to  ^-(n  irr  )  ,  i.e.,  the  change  in  particle  number  in  the 
necked-off  region.  This  subject  is  to  be  investigated  uner  1981  and 
1982  funding,  and  is  beyond  the  scope  of  available  theory  at  the  time 
of  this  report. 
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JUSTIFICATION  OF  DRIFT-SPEED-LIMITED  CURRENT 


When  the  ion  and  electron  temperatures  are  comparable  one  does 
not  get  appreciable  growth  of  the  ion  acoustic  instability  (Kovrizhnykh, 
1967;  Biskamp  and  Chodura,  1971),  but  the  modified  two-stream  instability 
(Ott  et  al,  1972;  McBride  et  al,  1972)  grows  with  linear  growth  rate 


Im  u>(MTS)  «-  t  wT11  »  i  u)  .  (1  +  u>2  /u.  2  )" 
2  LH  2  pi  pe  ce 


as  2  uiceV'A^Z 


where  is  the  "lower  hybrid"  frequency,  w  is  the  plasma  frequency, 

and  is  the  electron  gyrof requency .  This  persists  as  long  as  the 

drift  velocity  exceeds  approximately  the  sound  speed,  cg,  and  is 

marginally  stable  when  v  ~  c  .  The  nonlinear  saturation  time,  of 
^  u  09  S 

order  300  oj  ,  is  much  shorter  than  the  £/c  transit  time  of  electrons 
pe  s 

across  the  anode-cathode  gap,  so  the  waves  can  grow  to  saturation  within 
the  diode. 

Including  the  electron  magnetization  more  carefully,  one  also  gets 
an  "electron  cyclotron  drift  instability"  (Forslund,  Morse  and  Nielsen, 
1972)with  growth  rate 


Im  w(ECD)  r*  u>  — 
ce  v 

e 


where  v^  is  the  electron  drift  speed,  vg  is  the  electron  thermal  speed, 
and  the  electron  collision  frequency  as  before.  But  for  typical 
parameter  values,  v^/ve  should  be  of  order  7kfz  so  the  two  linear 
growth  rates  are  comparable: 


I.  u  -  2.6  X  10”  (fy  f)  \ 


I(MA) 
r  (mm) 


as  long  as  this  exceeds  in  Eq.  (1). 
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The  radial  gradients  in  magnetic  fields  can  also  give  rise  to  a 
"lower-hybrid  drift"  instability  (Krall  and  Liewer,  1971)  with  linear 
growth  rate 

„2  h 


j(LHD)  —  co  Nrr-  4^ - - —  (co  /co  ) 

v  1  pi  2B  dr  to  cl  ce'  pe' 
■*  ce  sJ 


7.4  x  10-5  co 


26  /10V 
A  \Z  ) 


hgB (MG) ‘ 


where  h. 


(ill) 


Since  B  is  at  most  2I/cr  ,  this  gives 

SL 


±  (4'2 


which  is  generally  less  than  the  MTS  and  ECD  instability  growth  rates. 


Because  the  modified  two-stream  instability  exponentiates  in  a 
few  picoseconds  and  saturates  in  several  growth  times,  it  is  most  likely 
that  the  electron  drift  speed  is  limited  in  this  way  to  values  near 
that  which  stabilizes  the  instability,  namely  the  sound  speed.  This 
is  the  justification  for  the  anomalous  drift-speed-limited  current- 
carrying  behavior  assumed  earlier. 
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APPLICATION  TO  CYLINDRICAL  IMPLOSIONS 
AND  THE  "SPLAT"  CODE 

From  Eqs.  (2)  and  (3),  we  easily  derive  a  fact  of  some  significance. 
Before  the  field  has  penetrated  a  cylindrical  z-pinch  plasma,  the  radial 
profile  of  the  Ez  electric  field  of  course  goes  from  negligible  values  at 
the  current  penetration  front  to  values  of  order  V/fi,  at  large  radii.  For 
all  radii  where 


» 


» 


I 


I 


E(r)  wv/”)  >  (f^)  (^(rf  (ir)(fo)5/2 


(33) 


the  current  is  drift-speed  limited.  For  temperatures  of  order  —  keV  and 
above,  the  critical  E(r)  value  at  r=a  is  much  less  than  V/i  'v  lMV/cm  when 


n^(r)  is  small  compared  with  1020. 


>-7 


But  from  the  voltage  rise  time  ^  10  sec,  we  can  estimate  a 
classical  skin  depth  in  the  dense  (classical)  region: 


Mc»>  -  29  \ev‘  (fr  fo)  - 


and  this  gives  a  core  current, 
a 

1^  *  J*  2"rrraE(a)e  ^a  r^^dr,  i.e. , 


(34) 


I1<MA)  -  0.21 


(  Tvt  /z_\2  t-1/4  (UnA)*  (26$  ni( 

yiQ-7 J  \10/  TkeV  \10  /  V A  /  a(nun)  101 


(a) 


(35) 


When  the  core  carries  most  of  the  current  (I»I^),  this  gives  for  the 
density  at  the  critical  surface  of  radius  a 


(36) 
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Thus,  for  a  keV  aluminum  plasma  carrying  1  MA  current  at  the  skin  of  a 
1  mm2  area,  if  the  on-axis  density  drops  below  about  1019  ions/cm3  any¬ 
where,  then  the  classical  "core"  tends  to  disappear  there.  And  since 
this  value  (Eq.  36)  of  n^  tends  to  be  small  compared  with  10  cm  ,  most 
of  the  E  field  is  shielded  from  the  classical  region,  i.e.,  E(a)<<V/Jl, 
while  there  is  a  classical  region. 

The  one-dimensional  strongly-coupled  plasma  implosion  code  SPLAT 
(Terry  and  Guillory,  1980b) ,  sets  (initially)  a  radius  a  beyond  which 
the  current  is  drift  speed  limited  at  the  initial  temperature.  In  the 
present  form  of  the  code,  this  radius  is  then  followed  in  a  Lagrangian 
manner,  moving  with  the  fluid  velocity.  In  fact,  however,  as  the 
temperature  and  current  increase,  the  real  "critical  surface"  at  which 
conduction  becomes  anamalous  moves  inward  more  rapidly  than  the  particles , 
following  the  variable  density  contour  of  Eq.  (36) .  Thus  a  larger  and 
increasing  fractionof  the  current,  Ij/I,  is  carried  "anomalously."  This 
fraction  is  thus  underestimated  by  the  code.  The  corona  is  actually 
more  important  than  indicated  by  the  code  for  another  reason:  its  density 
at  the  true  "critical  surface"  r=a  is  larger  than  presently  accounted 
for  and  the  total  ohmic  heating  should  be  larger. 
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SAUSAGE  INSTABILITY  DYNAMICS 


i 


I 


» 


» 


» 


» 


* 


The  nonlinear  behavior  of  the  axisymmetric  sausage  instability 
has  been  treated  only  in  a  few  limiting  cases.  A  nonradiating, 
perfectly  conducting  incompressible  MHD  model  with  uniform  density 
and  no  heat  conduction  was  treated  by  Book,  Ott  and  Lampe  (1976). 

In  both  the  long  wavelength  (1  >>  r)  and  short  wavelength  limits, 
these  assumptions  gave  tractable  equations,  which  showed  a  spool-like 
development  of  the  plasma  shape,  with  thin  rapidly  expanding  disks 
of  enlarged  radius  and  t  ,’oad ,  slowly  constricting  regions  of  reduced 
radius.  The  presence  of  finite  conductivity  and  axial  electric  field 
certainly  modifier  the  dynamics;  and  even  more  so,  the  inclusion  of 
finite  compressibility,  which  is  being  studied  now  by  Book  et  al,  may 
greatly  alter-*  the  exclusions  made  in  the  simpler  model. 

Because  of  chis,  we  have  proposed  in  our  10/81-10/82  program  an 
extension  of  the  Book  et  al  theory  to  include  finite  conductivity  and 
E  field,  and  the  construction  of  a  2-dimensional  r-z  code  for  modeling 
the  behavior  with  radiation,  circuit  model,  etc.  coupled  into  the 
dynamics. 


The  equations  describing  the  motion  are: 


continuity:  3tP  +  3Z(PVZ)  +  ~  3r(prvr)  -  0 


pDv/Dt  *  VP:  p(3  v  +  v3v  +  v  3  v  )  -  -  3  P 

t  z  zzz  r  r  z  z 


p(3  v  +  v3v  +  v3v)*=-3P 
tr  zzr  rrr  r 


state:  P  -  YTp/n^  +  Bz/8ir  +  KB2/r  B2  -  B2(r,z) 


(37) 


VxB: 

VxE: 


-3  (rB)  -  —  J  ;  -  3  B  -  —  J 

r  r  c  z  z  c  r 


3  E  -  3  E  ■  -  —  3  B 
z  r  r  z  c  t 


V*J»0: 


Fr  (rJr)  +  Vz  "  0 
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I  Heat  balance:  T  specified  as  function  of  discharge  voltage 

and  p,  in  simplest  model,  or  governed  by  time 
dependent  heating  and  cooling  rates,  in  more 
complex  model. 

I 

The  radial  component  of  J  cannot  be  neglected  because  then  9zB=0 
does  not  allow  the  instability.  The  constants  ka  in  the  unmagnetized 
conductivity  and  k^  in  the  Hall  conductivity  can  be  found  in  the 

I  standard  textbooks;  the  constant  kj  (see  Eq.  (4))  gives  the  drift-speed 

limited  current.  The  constant  K  provides  the  magnetic  hoop  stress  due 

to  B-line  curvature.  The  mass  density  is  denoted  by  o,  electron  (and 

ion)  temperature  by  T,  and  fluid  velocity  by  (v  ,v  ) .  All  field  and 

r  z 

*  fluid  quantities  are  functions  of  r  and  z.  Imposing  self-similar 

convex  profiles  of  P  with  radius  (e.g.,  Bennett  profiles)  and 
periodicity  in- z,  with  only  the  on-axis  density  and  Bennett  radius 
varying  in  z,  may  allow  certain  artificial  but  useful  simplifications. 

1  Making  the  ideal-MHD  assumption  is  not  admissible  here,  as  it  allows 

only  surface  currents  and  cannot  predict  the  density  drop  in  the 
necked-off  regions.  The  fluid  model  described  may  break  down  when 
these  regions  produce  non-stagnated  electron  flows. 


L 
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